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A search for supersymmetric partners of quarks is performed in the topology of multijet events 
accompanied by at least one tau lepton decaying hadronically and large missing transverse energy. 
Approximately lfb" 1 of pp collision data from the Fermilab Tevatron Collider at a center of mass 
energy of 1.96 TeV recorded by the DO detector is analyzed. Results are combined with the previously 
published DO inclusive search for squarks and gluinos. No evidence of physics beyond the standard 
model is found and lower limits on the squark mass up to 410 GeV are derived in the framework of 
minimal supergravity with tan/3 = 15, Ao = — 2mo and fj, < 0, in the region where decays to tau 
leptons dominate. Gaugino masses m 1 / 2 are excluded up to 172 GeV. 

PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 

Supersymmetric (SUSY) extensions [l[ of the standard squarks, and fermionic gluons (g), or gluinos, the super- 
model (SM) predict the existence of scalar quarks (q), or symmetric partners of quarks and gluons. These parti- 
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cles carry color and, if they are sufficiently light, they 
could be the most copiously produced SUSY particles 
at hadron colliders. The mass reach in direct searches 
at the CERN e+e~ Collider (LEP) is generally expected 
to be lower. However, direct searches for charginos and 
sleptons at LEP place stringent bounds in SUSY param- 
eter space which do translate into indirect constraints on 
squarks and gluinos [2]. At the Fermilab Tevatron Col- 
lider, previous searches for squarks and gluinos assuming 
R-parity conservation [J 0] were performed using events 
with jets accompanied by large missing transverse energy 
arising from the undetected lightest supersymmetric par- 
ticle (LSP) assumed to be the lightest neutralino (xT)- 
However, in some regions of the parameter space, squark 
and gluino cascade decays lead to final states with lep- 
tons. 

In the case of the third generation of sfermions, large 
mixing between the left and right handed chiral states 
may occur. Thus, in a given model, the lighter mass 
state of supersymmetric partners of the tau lepton (f-j^) 
might be the lightest of all sleptons [2]. In addition, if 
the is lighter than charginos (xt) & n d neutralinos 
(Xa)) they decay exclusively into staus, Xi ~~* t^v t and 
X2 — > t^t t , with the subsequent decay T t — > t ± Xi, 
leading to final states with tau leptons. Other decays to 
electrons and muons would be suppressed. In this region, 
sometimes refered to as the "tau corridor", squarks are 
lighter than gluinos and the production of squark pairs of 
the two first generations, pp — > QrQJj, dominates. Right 
handed squarks decay into a quark and the LSP while 
left handed squarks decay into a quark accompanied by 
the lightest chargino about 2/3 of the time or the sec- 
ond lightest neutralino about 1 /3 of the time. Figure Q] 
illustrates typical decay chains into final states with tau 
leptons. These cascade decays are characterized by the 
presence of two or more jets, large missing transverse en- 
ergy (JPt) and at least one tau lepton. Such signatures 
have not previously been exploited in SUSY searches at 
the Fermilab Tevatron Collider. 

In this Letter, we report on a search for final states 
containing at least two jets, at least one tau lepton de- 
caying hadronically and large missing transverse energy 
using 0.96 fb _1 of data collected in pp collisions at a cen- 
ter of mass energy of y/s = 1.96 TeV with the upgraded 
DO detector during Run II of the Fermilab Tevatron Col- 
lider. Squark production is investigated in the framework 
of minimal supergravity (mSUGRA) [5J with parameters 
enhancing final states with tau leptons. This analysis is 
further combined with the previously published DO anal- 
yses [4] developed in the framework of inclusive searches 
for squark and gluino production. 

The DO detector Q consists of a central tracking sys- 
tem, inside a calorimeter and a muon detector. The 
tracking apparatus is installed within a superconduct- 
ing solenoidal magnet of 2 T. It is composed of two 
subsystems, a silicon microstrip tracker and a central 




FIG. 1: Main Feynman graphs of left handed squark cascade 
decays into tau leptons; one tau lepton or two tau leptons are 
produced if squarks decay through a xt or a X2- 



fiber tracker, which provide efficient tracking and ver- 
texing over the pseudorapidity range \rj\ < 2.5, where 
T) = —In [tan [0/2)], 9 being the polar angle with respect 
to the proton beam direction z. The calorimeters use a 
sampling technology with liquid argon and depleted ura- 
nium organized in projective towers of size 0.1 x 0.1 in 
the T) — <f> plane, where (f> is the azimuthal angle. The 
central calorimeter (|r^| < 1.1) is separated from the end 
calorimeters which extend the pseudorapidity coverage 
up to ~ 4.2. Each calorimeter is located in its own cryo- 
stat, creating gaps in the interval 1.2 < \i]\ < 1.4, result- 
ing in a degradation in the energy resolution for jets in 
this region. The inter cryostat detector (ICD) is instru- 
mented with scintillating tiles. It provides additional cov- 
erage in the range 1.1 < \rj\ < 1.4. Calorimeter towers are 
subdivided into fine longitudinal layers. The ones closest 
to the interaction point form the electromagnetic (EM) 
section of the calorimeters, followed by the hadronic lay- 
ers which include the ICD. The muon detector surrounds 
the calorimeter. It provides coverage up to |?j| w 2.0. It 
is equipped with three layers of tracking detectors and 
scintillating trigger counters. An iron toroidal magnet of 
1.8 T is installed around the innermost layer. 

The data analyzed in this search were collected from 
April 2003 through February 2006. Events were selected 
online with unprescaled calorimeter based triggers de- 
signed to select dijet and multijet events with missing 
transverse energy. At trigger level, there are no require- 
ments on calorimeter clusters that would reject hadronic 
tau lepton decays, so these clusters will include both jets 
and tau candidates. The data sample corresponds to 
an integrated luminosity of (0.96 ± 0.06) fb _1 [7]. Due 
to the increasing instantaneous luminosity of the Fermi- 
lab Tevatron Collider, trigger requirements are tightened 
with time. In this Letter, we describe the most recent 
trigger requirements used. The DO triggering system 
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is organized in three successive levels (LI, L2 and L3). 
Events are selected by the LI if there are at least three 
trigger towers (of size 0.2 x 0.2 in 77 — plane) with trans- 
verse energy greater than 5 GeV. For dijet events, the 
modulus of the vector sum of the transverse momentum 
of simple cone jets, I/It, is required to be greater than 
20 GeV and 30 GeV at L2 and L3, respectively. The 
azimuthal angle between the two calorimeter clusters of 
highest transverse energy must be less than 168.75° and 
170° at L2 and L3, respectively. In addition, at L3, the 
scalar sum of the transverse momentum of the jets, Ht, 
must exceed 50 GeV and the minimal distance in azimuth 
between any calorimeter cluster and I/It is required to be 
greater than 25°. For multijet events, three calorime- 
ter clusters with transverse energy greater than 6 GeV 
and 20 GeV are required, and Ht must exceed 75 GeV 
and 125 GeV, at L2 and L3, respectively. At L3, I/It 
is required to be above 25 GeV. The trigger efficiency is 
greater than 94% for the signal events considered in this 
Letter that satisfy the offline selection criteria. 

We consider the signal event signature with two or 
more jets and at least one hadronically decaying tau lep- 
ton, accompanied by large missing transverse energy. SM 
processes involving W bosons decaying to a tau lepton 
and a neutrino, in association with jets, including top 
quark pair production and decay, contribute to the ir- 
reducible background. Another source of background is 
events with spurious I^t and electrons, muons or jets 
mimicking taus. SM expectations and signal efficiencies 
are computed with Monte Carlo methods, except for the 
multijet background from QCD processes, which is eval- 
uated directly from data. The detector geometry and re- 
sponse are simulated with a program based on GEANT 3 
[8j. One random beam crossing data event, recorded dur- 
ing data taking, is overlaid on each Monte Carlo event 
according to the instantaneous luminosity of the data 
sample in order to simulate detector noise and additional 
soft interactions. 

The hard scatter in SM processes is generated with 
version 2.05 of ALPGEN except for diboson and sin- 
gle top quark processes generated with PYTHIA 6.323 [loj 



and COMPHEP 4.1.10 [11(. The parton distribution func- 
tions (PDF) are modeled using the CTEQ6L1 H, Q]| li- 



brary. In all cases, the initial and final state radiation 
and the parton hadronization are simulated by PYTHIA. 
For all ALPGEN samples, the MLM parton-jet matching 
prescription [14( is applied to avoid duplicate phase space 
regions where both the matrix element as well as the par- 
ton shower contribute to the formation of parton jets. 
The inclusive tau lepton decay is simulated with the pro- 
gram TAUOLA 2.5 [15j with slighltly modified tau lepton 
branching ratios motivated by Next-to-leading or- 

der (NLO) QCD corrections for SM processes are taken 
into account by applying mcfm 5.1 17J NLO correction 
factors to the leading order (LO) cross sections calculated 
by the event generators. 



The analysis optimisation is performed in the 
mSUGRA model in a region of parameter space where 
final states with tau leptons dominate. While the uni- 
versal scalar and gaugino masses, mo and m l / 2 , are var- 
ied to explore the region of interest, a large mixing in 
the stau mass matrix is obtained by fixing the ratio of 
the neutral Higgs vacuum expectation values, tan f3, to 
15. The universal trilinear coupling, Aq, is set to — 2mo 
to favor high Higgs boson masses and the sign of the 
Higgs mixing mass parameter [i is negative so that the 
lightest chargino mass is increased while squark pair pro- 
duction cross sections are unchanged. This helps evading 
the indirect limits set by LEP. Squark and gluino pairs 
are produced, including all species but stop quarks, and 
their decay is simulated using PYTHIA. Sparticlc masses 
and couplings are computed with SUSPECT 2.3 [l8j], and 
SDECAY 1.1a [l9| . respectively. Samples are normalized 
to the NLO cross section computed with PROSPINO 2 (2pj 
and the CTEQ6.1M PDF set. The signal cross section is 
typically ~0.3 pb for squark masses of about 350 GeV. 

Collision data and simulated events are processed 
through the same reconstruction chain. Calorimeter tow- 
ers are clustered in jets with a cone algorithm of ra- 
dius TZ = v/(Ay) 2 + (A0) 2 = 0.5 [Hj], where y = 
^h\[(E + p z )/(E — p z )\ is the rapidity. Only jets with 
transverse momentum (pt) above 13 GeV are consid- 
ered. A jet energy scale correction (JES) is applied 
and simulated events are corrected to account for the 
measured jet reconstruction efficiency in data events. 
Hadronic tau lepton decays are characterized by a nar- 
row isolated calorimeter cluster with low track multiplic- 
ity 13, 23 1 . Tau finding starts with calorimeter clusters 
with Et > 5 GeV constructed with a simple cone al- 
gorithm with 1Z = 0.5. The innermost part of the cone 
(TZ = 0.3) is used to measure the tau candidate transverse 
energy while the energy deposited between TZ = 0.3 and 
TZ = 0.5 is used to compute the tau isolation variable. 
Tracks of pt > 1-5 GeV pointing to tau candidates in 
TZ = 0.5 are associated to them. In addition, we also use 
isolated tracks with pt > 5 GeV as seeds for taus. In 
this calorimeter cluster of TZ — 0.5 is constructed 

around the track and the rest of the tau construction 
procedure is the same as for the calorimeter seeded taus. 
All tau candidates are required to be associated with at 
least one track. One or two additional tracks, within 
2 cm in z of the highest pt track at closest approach, 
are added if the track invariant mass does not exceed 
1.1 GeV and 1.7 GeV respectively. Moreover, the third 
track is added only if the sum of the electric charge of the 
three tracks is equal to ±1. Electromagnetic sub-clusters 
with a minimum Et of 0.8 GeV are constructed from 
EM calorimeter cells belonging to the tau cluster. Tau 
candidates are further classified into three types, based 
on the tau signature in the detector. Type-1 taus have 
one associated track and no EM sub-clusters, typical of 
the decay into an isolated single pion or kaon. Type-2 
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taus have one associated track and at least one EM sub- 
cluster. Type-2 taus come mainly from decays via 
or af into a charged pion and one or more 7r°'s. These 
decay channels represent more than 50% of hadronic tau 
lcpton decays. They also include r — ► e decays. Type-3 
taus have more than one associated track, resulting from 
multiprong tau lepton decays. The tau candidate trans- 
verse energy is inferred from the transverse momenta of 
the associated tracks, the Et of the calorimeter cluster 
and the known calorimeter response to ir^ parametrized 
as a function of track px and rj. The IpT is computed 
from calorimeter cells prior to clustering. It is corrected 
for the JES and the pr of reconstructed muons. 

Preselected events are required to have at least two 
reconstructed jets and $t greater than 40 GeV. The 
event primary interaction vertex (PV) is restricted to be 
within 60 cm of the detector center along the z axis, 
I-Zpvj < 60 cm. The two jets with highest pt, ji and j2, 
are required to be reconstructed in the central part of the 
calorimeter (|?7| < 0.8), to have pr > 35 GeV and not to 
be collinear in azimuthal plane. Jet px can be mismea- 
sured if the primary vertex is incorrect, which can lead to 
apparent high IpT- This background from misreconstruc- 
tion, as well as instrumental backgrounds from noise in 
the calorimeter, is significantly reduced by requiring, for 
both jx and j 2 , CPF0 > 0.75, where CPF0 is the frac- 
tion of the sum of px of charged particles in a jet that 
point to the PV compared to the sum of pt of charged 
particles in the jet that point to any vertex. 

Events selected by the dijet trigger or the multijet trig- 
ger have different topologies. This fact is exploited by 
performing two different analyses, referred to as "tau- 
dijet" and "tau-multijet" , as summarized in Table [U 
In both selections, the IpT requirement is tightened to 
75 GeV. In the "tau-multijet" analysis, a third energetic 
jet is required and events with Ht less than 200 GeV 
are discarded. Multijet events with one mismeasured jet 
leading to IpT in the jet direction are rejected by select- 
ing events with IpT well separated in azimuth from the 
two leading jets. In the "tau-dijet" analysis, this proce- 
dure is extended to any jet reconstructed in the event. 
Major background contributions are W boson and top 
quark pair production, Z(—* ff)+jets events and multi- 
jet events. 

Events are required to contain at least one tau can- 
didate not already identified as one of the two highest 
Pt jets, ji and j2- Only tau candidates with Et above 
15 GeV are considered. The transverse momentum of the 
associated track in type-1 taus must exceed 4 GeV, and, 
in case of type-3 taus, the scalar sum of the pt of the asso- 
ciated tracks is required to be greater than 8 GeV. Quark 
or gluon jets are reconstructed as tau candidates with a 
probability between 0.2 and 0.8, strongly depending on 
energy and position in the detector. They are efficiently 
separated from hadronic tau lepton decays using Neural 
Networks (NNj), one for each tau type, which exploit the 
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for the "tau-dijet" and "tau- 


multijet" analyses. 




Preselection 


All analyses 


IpT 


> 40 GeV 


number of jets 


> 2 


A$(ji,j 2 ) 


< 165° 


[zpv| 


< 60 cm 


ji pt, M, CPF0 


> 35 GeV, < 0.8, > 0.75 


32 Pt, M, CPF0 


> 35 GeV, < 0.8, > 0.75 


Jets and IpT selection 


"tau-dijet" "tau-multijet" 


trigger 


dijet multijet 


IpT 


> 75 GeV > 75 GeV 


A$(#T,jl) 


> 90° > 90° 


A$(#T,j 2 ) 


> 50° > 50° 


A$ min (#T,any jet) 


> 40° 


number of jets 


> 3 


33 Pt, \rj\ 


> 35 GeV, < 2.5 


H T 


> 200 GeV 


Tau candidate selection 


"tau-dijet" "tau-multijet" 


number of tau candidates 


> 1 > 1 


A7?.(r can< j, ji) 


> 0.5 > 0.5 


A7?.(r can d, ji) 


> 0.5 > 0.5 


Optimisation 


"tau-dijet" OR "tau-multijet" 


S T 


> 325 GeV 


IpT 


> 175 GeV 



difference in transverse and longitudinal shower shape, as 
well as the isolation in the tracker and in the calorimeter. 
Training is performed with simulated Z — ► tt as signal 
and tau candidates in a multijet enriched sample from 
data as background. This analysis uses comparatively 
efficient NNj cut values [23j to ensure a high selection ef- 
ficiency for hadronic tau lepton decays (75%) while keep- 
ing a good background rejection (factor 14). The choice 
of moderate NNj requirements is in response to the very 
high rejection power of St and IpT selections against SM 



backgrounds, where St = Pt + Pr + -^f- Electrons are 
generally reconstructed as type-2 taus. Another Neu- 
ral Network was developed to separate electrons from 
hadronic tau lepton decays. In this case, Z — > ee events 
are used for background training. The cut applied gives 
a rejection factor on electrons of about about 20 while 
keeping 95% of hadronic tau lepton decays. Electrons 
in the ICD are reconstructed as type-1 taus. Therefore, 
type-1 taus in this region are discarded. Muons are de- 
tected as tracks in the central tracker and they leave 
a MlP-like signature in the calorimeter, which can be 
misidentified as type-1 taus. This contribution, as well 
as poorly reconstructed tracks, are suppressed for type-1 
and type-3 taus by comparing the energy of the calorime- 
ter cluster and the momentum of the associated track(s) 
(E/p > 0.7). In addition, the extremely high density of 
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the DO calorimeter favors the interaction of muons by 
Bremsstrahlung in the outermost layers. We reject tau 
candidates if they deposit more than 40% of their energy 
in those layers. This criterion also suppresses beam halo 
background and pions interacting by charge exchange in 
the calorimeter. 

At this stage of the analysis, events from W and top 
quark pair production account for 85% of the predicted 
background. According to the simulation, more than 
50% of tau candidates are due to a true hadronic tau 
at the particle level in all backgrounds. Non-simulated 
multijet background is also sizable in the "tau-multijet" 
analysis. It enters the sample when one jet mimics the 
signal tau signature. It is mainly distributed at low I^t 
as shown in Fig. [2] for both analyses. The background 
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FIG. 2: Distributions of JpT after the tau selection separately 
for the "tau-dijet" (a) and "tau-multijet" (b) analyses. The 
IpT cut is relaxed to 40 GeV. Data (points with error bars), 
SM predictions (plain histograms) with the estimated mul- 
tijet contribution added (white histogram) are shown with 
the (m ,m 1/2 )=(100,150) GeV signal point in the mSUGRA 
model with tan/3 = 15, Ao = — 2mo and /i < on top (dashed 
histogram). The multijet prediction is shown with ±1 stan- 
dard deviation (s.d.) error bands (light brown area). The 
uncertainties in the prediction of the other backgrounds are 
not shown explicitely. 



contribution from multijet events is estimated from data. 
For this purpose, the tau identification criteria on NNj 
for type-1, type-2 and type-3 taus are dropped to define 



a loose superset of data. The multijet contribution to the 
signal sample is estimated from the observed event yields 
in the loose and in the signal sample and from the rela- 
tive selection efficiency for events in the loose sample to 
also be part of this signal sample. The selection efficiency 
is known separately for the signal-like events, which we 
simulate, and for multijet background, which we do not 
simulate but estimate from data. The latter efficiency is 
estimated from a data control sample enriched in multijet 
events. To this end, only the $t requirement is changed 
to select low ]£t events (]pT < 75 GeV). This technique 
is also known as the "matrix method" . In the "tau-dijet" 
analysis, after the 1$t cut at 75 GeV is applied, the num- 
ber of events in each tau type is well predicted, as il- 
lustrated in Fig. [3J In addition to normalisation, the 
shapes of quantity are also well described by the sim- 
ulation. This includes observables used later on in the 
selection optimisation. For example, Fig. [3] shows that 
the Et of the tau candidate is well described by the sim- 
ulation, as well as St, & quantity indicative of the signal 
signature. Angular correlations are also well described. 
For example, the transverse mass of the tau candidate 
and IpT distributions are shown in Fig. [3] While the low 
edge of the distribution is shaped by kinematical cuts, 
it exhibits the Jacobian edge of the W boson transverse 
mass at higher values. This prevalence of tau candidates 
from W decay in background events, coming either from 
top quark pair or IT+jets production, is confirmed by 
the simulation. 

Finally, the "tau-dijet" and the "tau-multijet" selec- 
tions are combined with a logical OR, designated as the 
"tau" selection. The signal selection is optimized with 
two additional cuts on $t and on St, as quoted in Ta- 
ble |H The cut values are optimized by minimizing the 
expected upper limit on the cross section in the absence 
of signal computed with the modified frequentist CL a 
method 
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Statistical and systematic uncertainties are included 
in the CL S calculation. The systematic components are 
summarized in Table[TT] The JES uncertainty dominates: 
it is 25% for the background rate and 13% for the sig- 
nal rate. In comparison, uncertainties on jet resolution, 
on jet reconstruction and identification, and on the jet- 
vertex confirmation are negligible. Uncertainties related 
to tau leptons come from three sources and are small. 
The tau energy correction systematic uncertainty comes 
from the imperfect knowledge of the single pion response 
in the calorimeter in data. Measurements of the single 
pion response in data and their comparison with the sim- 
ulation result in an uncertainty of 6% in the single pion 
response used to derive the tau energy scale correction. 
This translates into a 1.2% uncertainty on the event yield. 
The NNj learning method is limited by the statistics of 
the data sample used for background. The corresponding 
systematic uncertainty is therefore estimated by fluctu- 
ating the NNj input variables according to the original 
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FIG. 3: Distributions of number of tau candidates per 
tau type (a), Et of tau candidates (b), St (c) and 
MT(j C and,IpT) (d) for events selected in the "tau-dijet" anal- 
ysis after the tau selection is applied. Data (points with 
error bars) and SM predictions (plain histograms) with- 
out the estimated multijet contribution are displayed with 
the (m ,mi /2 )=(100,150) GeV signal point in the mSUGRA 
model with tan j3 — 15, Ao = — 2mo and [i < on top (dashed 
histogram) . 



statistical uncertainty. The uncertainties from the NN 
against electrons are small in comparison. The uncer- 
tainty on the tau reconstruction and identification effi- 
ciency is dominated by the track finding efficiency. It is 
determined to be 3.0%. The systematic uncertainty on 
the integrated luminosity measurement is 6.1% Q and 
the trigger efficiency uncertainty on the signal and back- 
ground expectations is estimated to be 2% [4]. Based on 
MCFM, a 15% uncertainty is assigned to SM background 
cross sections. The uncertainty on the acceptance due to 
the approximate knowledge of the PDF is computed by 
using the CTEQ6.1M PDF error sets. The uncertainties 
in the modeling of the initial and final state radiation 
with pythia are estimated by varying the scale Q and 
the parton shower virtuality parameters [3] • 



TABLE II: Relative systematic uncertainties on SM back- 
ground expectations and signal events for the optimized "tau" 
analysis. 



Source 


Background (%) 


Signal (%) 


luminosity 


6.1 


6.1 


trigger 


2 


2 


jet energy scale 


25 


13 


jet resolution 


1.0 


1.0 


jet identification 


1 


1 


jet-vertex confirmation 


2 


2 


tau identification 


3.0 


3.0 


tau energy corrections 


1.2 


1.2 


tau AW selection 


1.2 


1.2 


cross section 


15 




PDF (acceptance) 


6 


6 


ISR/FSR 


6 


6 



Three data events are selected by the "tau" analy- 
sis. This is in good agreement with the SM expecta- 
tion of 2.3 ± 0.4(stat.) ± 0.7(syst.) events. Top quark 
pair production and W{— > ^)+jets events are the domi- 
nant background. The number of multijet events (-ZVqqq) 
selected by the "tau" analysis is estimated from data 
with the matrix method. Statistical and systematic un- 
certainties are included in the calculation, leading to 
Nqqjj = 0.1 ± 0.6 events. The number of events ex- 
pected for the signal point (m ,m 1 / 2 )=(100,150) GeV is 
6.5 ± 0.6(stat.) ± l.l(syst.). In this case, tau candidates 
are predicted to arise from hadronic tau lepton decays 
95% of the time. The distributions of $t and St are 
shown in Fig. [4] including the signal expectations assum- 
ing (mo,mi/ 2 )=(100,150) GeV. This point has a striking 
mass hierarchy where the f-j 4 is a few GeV lighter than 
the %2- I n this configuration, the lepton coming from the 
X° decay is mostly undetectable. Such mass configura- 
tions are difficult to detect in a search for direct xfx2 
production in the three lepton final states [26} . As shown 
below, this analysis is sensitive to this point at the 95% 
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C.L., whereas it is not excluded by the preliminary limits 
set on the slepton mass (mg > 100 GeV, mp, > 97 GeV 
and rrif > 93 GeV) and charging mass (m^± > 103 GeV) 
by direct searches at LEP2 [271 ]. 
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and factorisation scale are computed by varying by a 
factor two the default scale value Q in prospino. It 
is taken as Q — (rriq + mg)/2 and rrig, where rriq and 
rrig are the masses of the squark and the gluino, respec- 
tively for qq, qg and gg production. Renormalisation and 
factorisation scale uncertainties are typically 20% or less 
and PDF uncertainties can be as high as 40%. They are 
added in quadrature. The result is translated directly 
into a minimal observed limit and a maximal observed 
limit. They represent a band on the observed limit, com- 
puted with the default PDF and Q scale, as displayed in 
Fig. [5] In the case of the nominal signal cross sections, 
the "tau" analysis excludes squark masses, averaged over 
eight squark species, up to 340 GeV. The expected limit 
reaches the indirect limits from LEP2 chargino and slep- 
ton searches. 
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FIG. 4: Distributions of $ T (a) and St (b) with all the "tau" 
analysis requirements applied, except the one on the dis- 
played quantity. Signal events for the point (mo, mi/2)=(100, 
150) GeV in the mSUGRA model with tan/3 = 15, A = 
— 2mo and /i < are displayed on top of the background ex- 
pectations. The multijet background is found to be small and 
it is therefore not shown explicitly. 



From the number of data events, the SM expectation 
and the signal yields, upper limits on the production cross 
section are derived at the 95% C.L. The mSUGRA model 
with enhanced tau lepton final states is taken as a refer- 
ence and only the "tau corridor" is explored. Adding the 
estimated multijet event contribution to the SM back- 
ground expectation does not change significantly the re- 
sults. The multijet estimation is therefore not taken into 
account. Figure [5] displays the expected and observed 
limits in the mo— m 1 / 2 plane in the mSUGRA model with 
tan/3 = 15, Aq — — 2mo and \i < 0. In the region of inter- 
est, the analysis sensitivity is kinematically limited by the 
squark masses. Uncertainties on the signal cross section 
come from two major sources. The uncertainty from the 
PDFs is computed by summing quadratically the forty 
individual CTEQ6.1M error contributions. Uncertainties 
due to the imperfect knowledge of the renormalisation 
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FIG. 5: In the mo-mi/2 plane, expected and observed limits 
set by the "tau" analysis at the 95% C.L. in the mSUGRA 
model assuming tan /3 = 15, Ao = — 2mo and /i < 0. Limits 
are derived in the "tau corridor" only, which is represented 
on top of the regions excluded by LEP2 slepton and chargino 
searches. The observed limit is shown with a variation of the 
signal cross section by ±1 s.d. The minimal and maximal 
limits are the boundaries of the lighter area within the darker 
area of the "tau corridor" . 



The most recently published search for squarks and 
gluinos at DO |4j in events with jets and large miss- 
ing transverse energy using 2.1 fb -1 of data has not 
been specifically analyzed for tau lepton decays. How- 
ever, hadronic tau lepton decays can be reconstructed as 
jets and we have studied the sensitivity of these analyses 
to tau lepton final states and combined them with the 
present "tau" analysis to improve sensitivity. 

The published analyses are called "dijet" , "3-jets" and 
"gluinos" , each requiring an increasing number of jets. 
The "gluino" analysis requires four jets and it is not sen- 
sitive to the present signal topology. Regarding prese- 



10 



lection and, jets and Ipr selection, the "dijet" and "3- 
jets" analyses are identical to the "tau-dijet" and "tau- 
multijet" analyses, respectively (see Tablc|I|. The "dijet" 
and "3-jets" analyses veto high p? electrons and muons 
whereas the "tau" analysis rejects electrons and muons 
faking hadronic tau lepton decays. In the case of the "di- 
jet" and "3-jets" analyses, the optimisation is performed 
on the pair of requirements Ht and TpT- The optimized 
selections are H T > 325 GeV and $ T > 225 GeV for the 
"dijet" analysis and H T > 375 GeV and # T > 175 GeV 
for the "3-jets" analysis. Overlaps between selections 
are taken into account by defining a logical AND between 
them. While the "tau" analysis is limited to 0.96 fb _1 , 
the "dijet" and "3-jets" analyses include an additional 
1.17 fb _1 of data collected from June 2006 through Au- 
gust 2007. This leads us to consider seven exclusive chan- 
nels within the "early" dataset (0.96 fb _1 ) and three ex- 
clusive channels in the "late" dataset (1.17 fb _1 ). Table 
Mil details the ten channels of this combination and gives 
the number of selected and expected events for the SM 
backgrounds and for two signal points. Systematic uncer- 
tainties, as well as statistical uncertainties, are included 
in the CL S computation. 

Final limits are explored in the "tau corridor" in the 
mSUGRA model with tan/3 = 15, Aq — —2mo and fi < 
as shown in Fig. [6] The sensitivity of this combination 
reaches squark masses up to 408 GeV and the observed 
limit reaches 410 GeV for (m ,m 1/2 ) = (90,176) GeV or 
(mo,m 1 / 2 ) = (110,173) GeV. Although the "tau" analysis 
is performed on a dataset half the size of the "dijet" and 
"3-jets" analyses, Fig. [7] shows the relative gain in the 
production cross section upper limit is 10% if the "tau" 
analysis is included in the combination. This relative 
gain is 33% if datasets of equal integrated luminosity are 
considered by each analysis. 

The mass difference between the x 2 (or the xf ), the 
and the Xi vary inside the "tau corridor" , leading 
to different event kinematics. Model points close to the 
high mo border of the "tau corridor" predict that the 
f-j 4 mass is a few GeV below the % 2 mass, leaving the 
lepton produced in the x 2 decay with a low px- Going 
towards lower mo, slepton masses decrease and the mass 



difference xi ~ ^i increases, while the ff— Xi mass dif- 
ference decreases. Neighboring the LEP2 slepton limit, 
the point (mo,TO 1 / 2 )=(80,170) GeV exhibits similar mass 
differences between the x 2 an d the ff 1 (29 GeV) , and be- 
tween the fi and the xi (27 GeV). This point, with a 
squark mass of 396 GeV and a production cross section 
of 0.08 pb, is excluded by the "dijet", "3-jets" and "tau" 
analysis combination. We have explored these extreme 
mass configurations, as well as intermediate ones, and we 
exclude them. 

In conclusion, squark pair production was searched for 
in events with jets, tau lepton(s) and large missing trans- 
verse energy using 0.96 fb -1 of DO data recorded at a 
center of mass energy of 1.96 TeV during the Run II of 
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FIG. 6: In the mo-m^j plane, expected and observed lim- 
its set by the combination of the "tau", "dijet" and "3-jets" 
analyses. Limits are set at the 95% C.L. in the mSUGRA 
model with tan/3 = 15, Ao = — 2mo and fi < 0. Limits are 
derived in the "tau corridor" only, which is represented on 
top of the regions excluded by LEP2 slepton and chargino 
searches. The observed limit is shown with a variation of the 
signal cross section by ±1 s.d. The minimal and maximal 
limits are the boundaries of the lighter area within the darker 
area of the "tau corridor" . 



the Fermilab Tevatron Collider. A dedicated selection 
based on identified hadronic tau lepton decays has been 
developed. It was further combined with selections based 
only on jets and $t signatures and analyzing a superset 
of data that is twice as large (2.1 fb -1 ). No evidence for 
signal was observed in the combined analysis. The result 
of this search has been interpreted in terms of exclusion 
in the mSUGRA model with tan/3 = 15, Ao = -2mp 
and fj, < enhancing final states with tau leptons. The 
region of the parameter space where tau leptons are ex- 
plicitly produced, the so-called "tau corridor", is investi- 
gated. The highest excluded squark mass at 95% C.L. is 
410 GeV. 

This search is the first to explore supersymmetric mod- 
els in tau lepton final states in a multijet environment at 
the Fermilab Tevatron Collider. The combination of tau 
and inclusive analyses shows the highest sensitivity to 
the signature of interest. However, the "tau" analysis by 
itself is a crucial ingredient in order to obtain detailed in- 
sight into the nature of any signature from SUSY or other 
new physics that could soon be discovered at the Fermi- 
lab Tevatron Collider or at the CERN Large Hadron Col- 
lider. It is also important to complement other searches 
such as direct searches for chargino-neutralino produc- 
tion in the trilepton final state where one lepton could 
be kinematically undetectable. 

We thank the staffs at Fermilab and collaborating 
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TABLE III: Exclusive channels used in the "tau", "dijet" and "3-jets" analysis combination in the mSUGRA model assuming 
tan/3 = 15, Aq — — 2mo and /i < 0. "Early" and "late" refer to the two data taking periods used in this Letter. The flag "yes" 
or "no" means events passing the corresponding analysis are accepted or rejected for this channel respectively. The number of 
selected data, the SM expectations and the number of signal events for the mSUGRA parameters (mo, m 1( / 2 ) = ( 100,150) GeV 
and (mo,m 1 / 2 ) = (80,170) GeV are also given. The first quoted uncertainty is statistical and the second one is systematic. 



Chan. 


early 
"dijet" 


early 
"3-jets" 


early 
"tau" 


late 
"dijet" 


late 
"3-jets" 


Data 






SM 






Signal 
(100,150) 




Signal 
(80,170) 




1 


yes 


no 


no 






4 


4.2 


± 


0.5 


± 


1.1 


6.7 ± 0.6 ± 1.0 


3.5 


± 0.3 ± 


0.6 


2 


no 


yes 


no 






2 


3.9 


± 


0.4 


± 


0.8 


2.0 ± 0.4 ± 0.3 


1.5 


± 0.2 ± 


0.3 


3 


no 


no 


yes 






1 


1.1 
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FIG. 7: Observed and expected 95% C.L. upper limit on 
squark and gluino pair production cross sections with the 
combination of the "tau", "dijet" and "3-jets" analyses (la- 
belled "tau-combination" ) for mo = 80 GeV in the mSUGRA 
model with tan (3 = 15, Ao = — 2mo and fi < 0. The ex- 
pected cross section upper limit obtained with the "dijet" and 
"3-jets" analyses alone is also displayed ("jet-combination"). 
The NLO production cross section with a ±1 s.d. band, repre- 
senting the PDF, renormalisation scale and factorisation scale 
uncertainties is superimposed. 
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